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^ ! Abstract 
5h . 

CIh' Results of processing experimental data on charm production in hadron-hadron interac- 

tions are presented. The analysis is carried out within the frame of phenomenological model 
of diffraction production and quark statistics based on additive quark model (AQM). In low 
| energy region y/s = 20 — 40 GeV, the cross sections <j p n^cc.x (s), v-kN^ccX (s) are fitted by 

logarithmic function with the parameters connected by relationship of AQM. At collider en- 

t-H I ergies 200 , 540, 900, 1800 GeV, the values of a pp ^ c5 x{s) were obtained by a quark statistics 

method from the data on diffraction dissociation. It is established, that logarithmic function 

£^ | with universal numerical parameters describes the whole set of low-energy and high-energy 

data with high accuracy. The expected values of cross sections are a pp -, C cX = 250 ± 40 fib 
and 355 ± 57^6 at TEVATRON energy y/s = 1.96 TcV and LHC energy y/s = 14 TcV ac- 
cordingly. Opportunities of use of the obtained results for calibration of a flux of "prompt" 

o I ~ in high " energy component of cosmic rays are discussed - 

">< : 

1 Introduction 

The investigation of charm production processes at overaccelerator energies is an actual 
problem of elementary particles and cosmic rays physics. Nowadays there are four groups 
.£h ! of experimental data: 

1) data on total cross sections of charm production in hadron-hadron interactions 
c3 ' vtotipN — ► cc + X), <j tot (irN — > cc + X) in the region of low energies y/s = 20 — 40 GeV; 

2) data on total cross sections of charm production in a photon - hadron interactions 
VtotilN — > cc + X) at low energies y/s = 6 — 20 GeV and - at virtualities of a photon 
up to Q 2 = 4 GeV 2 - in the region of high energies y/s = 160 — 240 GeV; 

3) data on production of J/^f, *$f(2S) particles in pp-interactions at y/s = 1800 GeV 

4) data on differential cross section of charm production in pp— interaction in the 
central area at y/s = 1960 GeV and transverse momentum pt > 5.5 GeV 

The first two groups of the data and the appropriate review of the literature contain, 
for example, in Ref. jlj. Only the first group of the data has the status of a direct 
information on charm production in hadron-hadron interactions Three other groups differ 
from the first one by initial or final states, therefore the joint model independent processing 
of all three groups of the data is impossible. Accordingly, comparison of the theory to 
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experiment is at a loss also: fixing of parameters of theoretical model by its comparison to 
the one group of high-energy data results in catastrophic disagreements of the same model 
with other group of the data. In the scientific literature this situation is characterized as 
"an unsolved problem" 

In the present work an estimation of the total charm production cross section is offered. 
The estimation is based on processing of collider data on diffraction dissociation with use 
of phenomenological model of diffraction charm production and quark statistics. The 
result reduces to simple formula with well enough certain numerical parameters: 

g 

&pp/pp—tcc+x(s) Cpp/pp In ) 

so (1) 
C p p/ pp = 26.78 ± 1.44 //&, ^ = 18.23 ± 0.33 GeV. 

For TEVATRON energy ^ = 1800 GeV our estimation gives a Pp ^ 5c+x ^ 246 ± 30 pb; 
for LHC energy ^Ts = 14000 GeV, (r pp ^ c5+ x =^ 356 ± 35 fib. 

In cosmic ray (CR) physics charmed particles, arising in interaction of primary particle 
with an atmosphere, are sources of so called "prompt" high-energy muons, carrying away 
energy from Extensive Air Shower (EAS). A well known problem about the nature of a 
break ("the knee") in the cosmic rays energy spectrum at Ek(iab) ~ 3 x 10 3 TeV {^fs~k ~ 
2.5 TeV) can not be solved without the quantitative information on the contribution 
of charmed particles in the total muon flux. Obtaining such information is a complex 
problem which is solved by joint processing the accelerating data (LHC data will be 
added to them in the coming years) and results which should be received at CR facilities. 
These are accelerating data on total and inclusive cross sections of hadronic interactions 
(including cross sections of charm production) and CR data on the flux of muons with 
energies E^ Lab) > 70 TeV. 

In conditions of incompleteness of accelerating data, the measurements of "prompt" 
muons fluxes at energies E p > 70 TeV gain special interest - their results are capable to 
give the preliminary information about a character of hadronic interactions in the range 
of ultrahigh energies. For planning and preparation of the experiment, the forecast of 
expected value of charm production cross section is necessary, provided that laws of this 
process are extrapolated from accelerating energy region. Our estimation represents one 
of the variants of such forecast which can be used for calibration of the flux of ultrahigh 
energy CR muons. 

2 Diffraction model of charm production and additive quark 
model 

The offered method of the estimation of charm production cross section is based on the us- 
age of additional set of experimental data. The question is cross sections of proton diffrac- 
tion dissociation in pp - interactions at SPS CERN energies a/s = 200 GeV, 900 GeV 
and TEVATRON ones yfs = 546 GeV, 1800 GeV 7J. According to these data, the 
diffraction dissociation cross section depends logarithmically on the energy: 

a DD {pp^X) = C DD \n-. (2) 

so 

The basic assumptions consist in the following: 
1) charm production occurs in process of diffraction dissociation; 
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2) charm production cross section is extracted from total cross section (J2J) by quark 
statistics rules obtained from additive quark model: 



(T t ot(pP ^cc + X)~ k c5 x a DD (pp -> X), k c5 ~ 0.025 ± 0.004. (3) 

These assumptions, from our point of view, have general enough phenomenological 
character. Formulas (J2J), © have the status of the asymptotic ones. They are fulfilled 
more precisely, the greater is the energy of interaction. The value of k c5 — 0.025, presented 
in (jHJ), is obtained from the relationships in Ref. [Bj: 

uu : dd : ss : cc = 1 : 1 : 0.38 ± 0.07 : 0.06 ± 0.01, (4) 

here the symbol qq designates the number of secondary hadrons containing quark q or an- 
tiquark q. Existence of similar relationships is predicted by additive quark model (AQM) 
jU]. Direct check of these relationships is realized in experiments, in which the composi- 
tion of secondary particles is supervised. A condition of inclusion of charmed particles in 
quark statistics is an absence of energy restrictions on process of their production. This 
condition is obviously fulfilled for charmed particles produced in the diffraction region. 
All numbers appearing in (j3J), except the one concerning to c— quarks, are obtained by 
comparison of AQM to the experiment at y/s = 540 GeV ^U]. An estimation of the sup- 
pression factor of charmed hadrons A c /A s ~ 0.15 is based on characteristic assumption 
for AQM: for q = s, c, b 

% l 

where the summation is over all particles and resonances with spins J and masses M, 
containing a quark q or an antiquark q. Each of the statements reflected in formulas 
(J2J) - © was repeatedly discussed in scientific literature. We shall note, that typical 
quantitative accuracy of the phenomenological models based on these statements is 10 - 
20 %. 

First of all we shall find out the opportunities of AQM and logarithmic dependence for 
cross sections in the description of direct experimental data on charm production in pN 
and 7iN interactions at energies = 20 — 40 GeV. These two hypotheses are reduced 
to formulas 

s 2 

&pN-+cc+x(s) — Cpjvlri — , (TTrN^cc+x{s) = ~ a pN-+cc+x(3 s /2), (5) 

So O 

containing only two free parameters. Numerical values of these parameters are defined by 
joint fitting of pN and data: 

C pN = 28.84 ±2.10 fjb, ^ = 18.51 ± 0.36 GeV, x 2 = 0.89. (6) 

Results of fitting are shown in Fig. 1 (ttN data are shifted up on 20 fib for clearness). 
As it follows from Fig. 1, the phenomenological diffraction model of charm production 
describes well the direct low-energy data. However extrapolation of these results in the 
region of high energies ^/s = 1000 — 3000 GeV has no substantiation, certainly. 

As it was already marked in the beginning of the Section, there is an opportunity 
to involve the additional experimental information in the range of SPS - TEVATRON 
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Figure 1: Accelerator data on charm production cross section in pN and ttN interactions. ttN 
data are shifted up for clearness. 



energies. Cross sections of diffraction dissociation measured at these colliders have the 
following values: 

y/s, GeV 200 546 900 1800 



a DD} mb 4.8 ±0.5 7.89 ± 0.33 7.8 ± 0.5 9.46 ± 0.44 

It is easy to verify that the given numbers are fitted in logarithmic dependence. Cal- 
culation of values of charm production cross sections from the formula Q, taking into 
account relationships of AQM (HJ), gives: 

y/s, GeV 200 546 900 1800 



a pP ^ c5+x , lib 120 ± 21 197 ± 32 195 ± 32 236 ± 38 

The obtained values of charm production cross sections have the status of model 
dependent processing collider data. 

In Fig. 2 these results are presented together with the data of direct measurements in 
the region of low energies. The result of fitting the total data set is: 

C pP = 26.78 ± 1.44 /jjb, yfa = 18.23 ± 0.33 GeV, x 2 = 0-98. (7) 

As it is evident from the foregoing, values of the parameters © and (J7J) coincide within 
the limits of statistical errors. It means the charm production cross section in hadronic 
interactions is described by universal logarithmic dependence covering both low-energy 
and high-energy experimental data. 
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Figure 2: Joint fit of low-energy and high-energy data of charm production. 



3 Discussion of results and conclusions 

Only the accuracy level is surprising with which low-energy (JHJ) and total (J7|) parameters 
of logarithmic function coincide. With regard to ideas of phenomenological models of 
diffraction production and quark statistics, they were repeatedly tested at processing 
most various experimental data. These models give reliable quantitative predictions at the 
accuracy level 10 - 20 % after fixing several parameters by experimental data. Therefore 
we consider possible to use the results (|7j) for forecasting cross section of charm production 
in pp/pp interactions in overaccelerating energy range. 

In Fig. 3 the dependence of charm production cross section on energy is presented. In 
Figure is shown the dependence with parameters (J7|) extrapolated in overaccelerating 
energy range. The confidence strip of the forecast corresponds to 90 % CL. The expected 
value of cross section is a pp ^ C c + x — 270 — 300 fib in the energy region ^fs ~ 3000 — 
5000 GeV and a pp ^ c - c+x ^ 355 ± 35 fib at LHC energy ^1 = 14000 GeV. 

Using the quark statistics, the total cross section of charm production can be presented 
as the sum of inclusive cross sections. Appropriate relationships between cross sections 
^] are checked experimentally up in the region y'i = 20 — 40 GeV: 

a(D s D s ) /a(D°D° + D+D~) ~ 0.2, a(A c D) /a(D°D° + D + D~) ~ 0.3 

As it was shown above, a hypothesis of scale invariance of total cross section parameters 
is experimentally motivated; therefore we expand this hypothesis on the parameters of 
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Figure 3: Extrapolation of the charm production accelerator data fit to overaccelerating energy 
range. 

inclusive cross sections. In this case we have in laboratory system 

a DD (E) « 154.1 + 17.86 In E fib, 

°d 3 d s (E) « 30.8 + 3.55 In £ fib, (8) 

a AcD (E) ^ 46.2 + 5.33 In E fib. 

where E is in PeV (1 PeV= 10 3 TeV). A technique of use of similar type formulas for 
calculation of an expected flux of "prompt" high-energy muons, accompanying EAS, is 
described in Ref. jllj . Detection of such muons, in essence, is check of numerical values 
of factors appearing in formulas (jHJ). 

Now the experiment on recording very high energy (VHE) muons (E^ > 70 TeV) 
is carried out at Baksan Underground Scintillation Telescope (BUST) [T2|. Preliminary 
results of this experiment, expected in the near future, will show as far as the offered 
estimation of charm production cross section is adequate. It should be noted that limiting 
energy of TEVATRON y/s max = 1.96 TeV is very close to the energy of the break in CR 
energy spectrum ^Sk nee ^2.5 TeV. If total cross sections of charm production would be 
measured at TEVATRON, their small extrapolation would allow to obtain the information 
necessary for calibration of the flux of VHE muons. In this case, logic and interpretation 
of the experiment on recording CR muons would get at once a rigid character. 
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It should be pointed out, there are three possible (and alternative) variants of results, 
each of which represents the certain interest. 

1) The measured flux of CR muons within the limits of measurement errors coincides 
with calibration flux calculated according to the estimation of charm production cross 
section (0). In this case the formula (JTJ can be used for testing microscopic models. 

2) The measured flux of CR muons will exceed noticeably the calibration value, but 
the appropriate cross sections of charm production at energies in the region of a break of 
CR energy spectrum will stay essentially smaller 10 mb. In this case it will be necessary 
to recognize, that either quark statistics rules, in that kind in which they were used in 
the present work, become incorrect in the area of high energies, or the new mechanism of 
charm production, distinct from the diffraction one, is included at these energies. 

3) The measured flux of CR muons will correspond formally to charm production cross 
section exceeding 10 mb. This result will mean that VHE muons carry away the energy 
from EAS forming an observable break of CR spectrum (see, for example, [T3*]). However, 
it is necessary to note, that sources of such big numbers of VHE muons, most likely, are 
not reduced to charmed particles - a New physics will required for interpretation of such 
effect. Certainly for the formulation of new ideas, the data on the flux of CR muons 
should be considered together with all other accelerating and CR data. 
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